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Abstract. Recently efforts have been made to use generalized sine functions to perfectly 
reconstruct various kinds of non-bandlimited signals. As a consequence, perfect recon- 
struction sampling formulas have been established using such generalized sine functions. 
One of our goals in this paper is to study the error when a truncation scheme is employed 
for this type of sampling formula which involves an infinite sum. This will be done in Sec- 
tion [3l We further analyze the error of the sampling formula when there are noises present 
in the samples in Section 2J Finally we discuss the reproducing properties and Sobolev 
smoothness of functions in the space of non-bandlimited signals that admits such a sam- 
pling formula. We begin with a discussion of properties of the generalized sine function in 
Section [2 



We begin by establishing some notations used in the paper. Let N be the set of natural 
numbers, Z be the set of integers, and Z + := {0} U N. For a positive integer n 6 N, we use 
the index set Z n := {0, 1, . . . , n — 1}. Furthermore, we denote by R the set of real numbers, 
and by C the set of complex numbers. Let X be a subset of R, and for q 6 N, we say a 
function / is in L q (X) if and only if 
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1. Introduction 




and / is said to be in 



L°°(X) if 



11/IU.x := esssup{|/(t)| : t G X} < oo. 



l 



Similarly, let Z be a subset of Z, a sequence y '■= (yk '■ k e Z) is said to be in if and 

only if 

I \ ^ 

\\y\\i,z '■= I z~2 I < 00 > 

and y is in l°°(Z) if 

||?/||oo,z := sup{|y fc | : k G Z} < oo. 
In digital signal processing, the classic sine function is fundamentally significant due to 
the Whittaker-Kotelnikov-Shannon (WKS) sampling theorem [161 El 12] • Recall that the 
classic sine function is defined at t G R by the equation 

, . sint 
sine (t J := . 

The WKS sampling theorem enables to reconstruct a bandlimited signal from shifts of 
sine functions weighted by the uniformly spaced samples of that signal. It is natural to ask 
whether similar sampling theorem exists for non-bandlimited signals. To that end, recently 
efforts have been made to extend the classic sine to generalized sine functions, for example, 
in |H El E]. One kind of generalized sine functions given in [I], denoted by sinc^, is defined 
as the inverse Fourier transform of a so-called symmetric cascade filter, denoted by H. 
Specifically, 



(1.1) sinc F := d^' l H, 

where for any sig nal / G L 2 (R) and £ G R 

= /(£):= 4= / f(t)e-^dt. 



The symmetric cascade filter H is a piecewise constant function whose value at £ G R is 
given by 



;i-2) H(0 := bnXiAO, 



where the sequence b = (b n : n G Z + ) is in / 2 (Z + ), ^/ is the indicator function of the set /, 
and the interval I n , n G Z + , is the union of two symmetric intervals given by the equation 

I n := (-(n + l),-n]U[n,(n + l)). 

Of course, we have that if G L 2 (M) because 6 G / 2 (Z + ), and hence sinc^p G L 2 (R) since the 
Fourier operator is closed in L 2 (1R). 

With the generalized function sinc^, a perfect reconstruction sampling theorem was 
established in [3] for the purpose of reconstructing non-bandlimited signals. This kind 
of reconstruction sampling theorem may be very useful to study signals with polynomial 
decaying Fourier spectra that arise in evolution equations and control theories [131 IT3] . 

One of our goals in this paper is to study the error when a truncation scheme is employed 
for the above sampling formula which involves an infinite sum. This will be done in Section 
|3j We further analyze the error of the sampling formula when there are noises present in the 
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samples in Section [H Finally we discuss the reproducing properties and Sobolev smoothness 
of functions in the space of non-bandlimited signals that admits such a sampling formula. 
We begin with a discussion of properties of the function sinc^ in Section [2j 



2. Properties of the generalized sinc functions 

Surprisingly the function sincj? has many properties that are similar to the classic sinc, 
such as cardinal, orthogonal properties and it behaves also similarly to the classic sinc. 
In the special case, the function sinc^ reduces to the classic sinc. Let us first review the 
approach to obtain an explicit form of sinc^. 

The symmetric cascade filter H can be associated with an analytic function F on the 
open unit disk 

A:={CGC:|C|<1}. 
The value of F at z G A is well defined by 

(2.1) F(z) := 

nez+ 

as b G / 2 (Z + ). Recall that the Hardy space H 2 (A) consists of all functions / analytic in 
A, with norm given by 



re(0,l) ^ Jf-Tr.vrl 



re(0,l) Z7r </[-7r,7r 

Since, by hypothesis, b G Z 2 (Z + ), we have that F G H 2 (A). Consequently its extension to 
the boundary dA of A is in L 2 (dA). 

Thus, from equations (11.11) . (II. 2p and (12. ip an explicit form of sincp^), t G R can be 
found as 



(2.2) sincF(t) = sinc Qj Re ji^e^e^*} , a.e. 



where Ke(z) is the real part of a complex number z. 

We observe that if b G then H G L 1 (M) and F is continuous on the boundary of 

A, which in turn implies sinc^ is continuous and bounded. 

A very interesting fact, as discovered in the paper [3], is that when F is imposed with 
a stronger condition of having analyticity in a neighborhood of the closed unit disc A, the 
function sinc^ can be generated through a function, denoted by G, that is also analytic in a 
neighborhood of the closed unit disc, real on the real axis and normalized so that G(l) = 1 
and G"(l) ^ 0. The function G is linked to F by the equation 

(2.3) m == 

A real- valued function (fie whose value at t G R is then defined through the imaginary part 
of the values of G on the unit circle by the equation 

(2.4) Mt ) := 



Applying equation (12. 2p to compute sinci? by using equations (12.31) and (12.41) yields for 
all t £ E, 

(2.5) sinc F (t) = (j) G {t). 

Two important examples can be demonstrated for this construction. When G = z, i.e., 
F = 1, we have sinc^ = 4>g = sine. For the second example, let G be the Blaschke product 
of order tiGN with parameters a := (aj : j £ Z n ) £ (—1, l) n , that is, 

GW = a. w := n 



Then sincp(£) = (f)c(t) = sm9 "^ ; where # a is determined by the boundary value of the 
Blaschke product at i £ K by e Wa ^ = B a (e u ). 

We next list some properties of the function sinc^. 

Proposition 2.1. Let the generalized function sinci? be defined by equation (12. 2p . Then 

(1) sincF(ri7r) = F(l)5„ i0 , where b~ n fl = 1 if n = and <5 nj o = if n E Z \ {0}. 

(2) sinc^ is bounded, infinitely differentiate. 

(3) |sinc F (i)| < 4 " b ^ 1 | ^ +) , for t £ R, and sinc F £ L 2 (R). 

(4) TTie set {sinci?(- — nit) : n £ Z} is an orthogonal set, that is 

(sincF,sincF(- - rwr)) = 7r|| b||p( Z+ )^o, 

where (•, •) is t/ie usim/ inner product on the Hilbert space L 2 (R). 

Proof. The first two statements directly follow from equation (12.21) . The third statement 
follows from equation (12. 2p and noticing sinc(t) < for any £ £ R. The fourth statement 

is a special case of Corollary 3.2 of [4J. For the convenience of readers, we provide a direct 
proof here. By Parseval's theorem and equation (II -ip we have 

I sinc F (t)sinc F (t-mr)dt = - / # 2 (a;)e in7ra: ds = J / J" b 2 kX i k (x)e in7rx dx 
Jr 2 J k 2 J M k&+ 

( \ 

where, in the last equality we have used the orthogonality of the set {e _m7r ^ : n £ Z} on 
Jfc, k £ Z + . The interchange of the integral operator and the infinite sum is guaranteed by 
the convergence of the series. 

□ 

3. Sampling Truncation Error Analysis 

In [I], a Shannon-type sampling theorem is given concerning functions in the shift- 
invariant space 



e mnx dx 



TV 



Inez 



\ /, : { > v C n sinc F (- - nir) : F(l) = 1, c = (c n : n £ Z) £ Z 2 (Z) 
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of sinc^. The Shannon-type sampling theorem is the direct result of the properties in the 
previous proposition. We record it here. 



Theorem 3.1. A signal f G Vp if and only if 

(3.1) / = ^/(rar)sinc F (- - mr). 

Equation (13.11) necessarily implies that the sampling sequence (f(n7f) : n G Z) G / 2 (Z) 
by the orthogonality of the set {sinc F (- — nn) : n G Z}. The above equation of course is 
true in L 2 (R) norm. However, if b G l l (7* + ), equation (13.11) holds true pointwise, because 
by Cauchy-Schwartz inequality the series on the right side of equation (13.11) converges 
uniformly, hence the limiting function / is continuous. 

We remark that, as pointed out in [3], a function / G Vp can be characterized by its 
spectrum. Specifically, a function / G Vp if and only if 




(3-2) 77(£) = \J — ( ^/(mr)e- 1 ^ H(£). 

Equation (Q holds true in L 2 (R) if b G l 2 (R), and a.e. pointwise if / G L X (R) and the 
sample sequence (f(nir) : n G Z) G / X (Z). 

The following property is true for functions in the space Vp that is similar to Parseval's 
identity. 

Proposition 3.2. If f E Vp then 

(3-3) ll/lli 2{K ) = ^ll&ll^ + )E/ 2 (^) 

Proof. This is a direct result of equation (I3.ll) and Proposition 12. II (4) . □ 

The recovering formula in ( 13. lft involves an infinite sum. In practice, we need to truncate 
the series to approximate /. Here, we prefer an adaptively truncated sum and offer a 
pointwise estimation of the error. 

For fixed n G N and t G R, define the index set 

: = {j -jeZ,\t- jn\ < nir} 
and the partial sum given at t G R by 

(3.4) S n (t):= /(jvr)sinc F (t-jvr). 

The adaptive truncation strategy allows that for any t G R and n G N, there are approx- 
imately 2n functions sinc F shifted by a distance of jtt, j G J n {t), on both sides of f. The 
following theorem states that the truncated error estimate is 0(n -1 / 2 ). We recall a Calcu- 
lus fact that will be used a couple of times later. For a positive and decreasing sequence 
a := (a n : n G Z + ), if ^— j- — ^- = c, where c is a positive constant, then 

/ \ 1 11/" 1 



" . a o C J [ao.oo) 37 



Theorem 3.3. Let f G V F and b G / 1 (Z + ). Then we have 

\b 



(3.6) 



sup\f(t)-S n (t)\ < 



,|&||l 2 (Z+) 

Proof. We first note that by equation (I2.2p . for t G 



-T 1 -4+- 1 



7T 3 \ n 2 n 



(3.7) 



|sinci?(t)| < ||6||zi 



sinc- 
2 



< 



1% 



|t/2| 



if lil < 2, 



if lil > 2. 



By Cauchy-Schwartz inequality and equation (13.71) . we have 



\f(t)-s n (t)f 



f U^smc F (t - 

j€Z\J n {t) 



< 



E ^ 

v i6Z\J„(t) 




sinc|(t-i7r) 

j<EZ\J n (t) 



(3.8) 



< 4||&||, 2 1(z+) 




E 

jez\J n (t) 



(t - jTif 



We next estimate for t G M the value of 
(3.9) R,,(t) := 



E 



1 



(t-jn) 2 ' 

It is easy to see that i?„ is periodic with period 7r. For £ G [0,7r), using equation (13.51) we 
have 



(3.10) 



1 



1 



Rn(t) < 2 ( , 

Using equation (13. 3p . we obtain that the quantity 
(3.11 



-4dx | = — r ( — r + — 



E r-w<T.f 

jez\J n (t) 



n-K) 



7r 2 \ n 2 n 



2 

i 2 ( 



Finally we conclude (13. 6 \ by combining equations f 13 . 8 1) . (13 . 1 0[) and (13. lip . 

4. Error analysis when noises present 



□ 



Theorem 13.11 establishes that a signal in the space Vp can be perfectly reconstructed 
by an infinite sum of shifts of the generalized sine functions weighted by equally spaced 
samples of that signal. However, samples are often corrupted by noise in practice. In 
[15] . Smale and Zhou gave an error estimate in the probability sense for Shannon sampling 
theorem with noised samples. In [I], Aldroubi, Leonetti and Sun studied the error by frame 
reconstruction from noised samples. In this section, we shall investigate the error of the 

6 



sampling formula (13. ip with noised samples. Specifically, we deal with the following noise 
model concerning the noisy samples f(nir), nfZ whose value corrupted by noise is given 
by 

(4.1) /(mr) = /(mr) + e(rwr), n e Z, 

where we assume that (e(mr) : n G Z) is a sequence of independent and identically dis- 
tributed random variables with the expectation and variance of each given by 

(4.2) E(e(mr)) = 0, Var(e(mr)) =a 2 , n G Z. 
Thus in practice, we recover / G Vf by 

(4.3) /■> = 7(n7r)sinc i r(- - mr). 

We next study the expectation E(/ — p) and variance Var(/ — Z' 1 ). 

Theorem 4.1. Let f & Vf be recovered by equation (14.31) wit/i t/ie noised samples f(mr) 
being referred to in (14. ip . T/ien 

E(/(t) -/"(«)) = 

and 

(4-4) Var(/(t) - /1(f)) < 2 ( r 2 ||6||f 1(z+) (l + . 

Proof. We first compute the expectation E(/(t) — f\t)). 

E(/(t) - /*(*)) = E(]T e (mr)sinc F (t-mr)] 

= E(e(mr))sincF(f — mr) 
= 0. 

Invoking the assumed independence of e(jir), j G Z we obtain 

Var(/(t) - f\t)) = Var ( ^ e(mr)sinc F (f - mr) J 

VngZ / 

= Var(e(n7r)sincj7(t — mr)) 



nG2 



(4.5) = a 2 y^sinc|i(£ - mr). 
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The sequence sinc^(£ — mr) can be easily estimated. Note it is periodic with period ir. 

Recalling equation (13. 7ft . for t G [0, tt), we have 

sinc^(i — mr) = sinc^(t) + sinc^(i — tt) + sinc^(£ — ri7r) + sinc^(i — n7r) 

neZ ne-N nel+N 

< 2+ sinc^(£ — mr) + sinc^(t — n-yr). 

ne-N nel+N 

Noting for £ G [0, 7r), n G N, ^~ 2 W7r ^ > | and in view of equation (I3.5P we obtain that 



2/ n .1,1.2 ( 4 2 f 1 

^ 8 inc F (t-n7r)<||6|| I1(z+) I - + -/-, 

nel+N \ -M^' 00 ) 



3 dx ) = H b ll^ + ) ( ^2 ) " 



Similarly we have 



ne-N 

Consequently we obtain that 



smc 2 F (t-mr) < \\b\\ 2 1{z+) (^j . 



(4.6) 



^sinc|(t-mr) < 2\\b\\l (z+) M + — J . 



Finally combining equations (14. 5 j) and (14. 6 p proves equation (14. 4p . □ 

5. The reproducing property and Sobolev smoothness 

When the analytical function F is chosen to be F = 1, the space Vp reduces to the space 
of bandlimited signals. The space of bandlimited signals is a reproducing kernel Hilbert 
space (r.k.H.s) [TJ. We next show that the space Vp has a similar property. However, the 
reproducing kernel is a distribution in the space of tempered distributions. We define the 
distribution for x, t G R by 

$(x,t) := -smcp(t - x) ^ e i2kx . 

Recalling the Poisson formula in the distribution form: 

^ e -i2fc« = ^J^Six - kir), 
fcez fcez 

where 5 is the usual Dirac delta function, we immediately obtain an alternative form of $ 
given by 

(5.1) $(x, t) — 5{x - kn)smc F (t - x). 

fcez 

Theorem 5.1. Let f G Vf then 



(5.2) /(*)= / f(x)$(x,t)dx 

Jr 

in the distribution sense. 
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Proof. Poisson's summation formula indicates that 



(5.3) £/(mr)e 



fcez 




This leads to 



/(*) = r- 1 [J2 j: f(-+M) H (*) 



= ^^- 1 (^/(- + 2A;)i7)(t) 



The interchange of the order of the sum and integral operator is justified by the convergence 
of the series. Recalling equations fll.ip and the convolution theorem for Fourier transform, 
we therefore conclude that 



Next we discuss the Sobolev smoothness of a function in the space Vp. We say that a 
function / belongs to the Sobolev space H S (M.) if 



The Sobolev smoothness of / is defined to be v%{f) '■— sup{s : / 6 iiP(R)}. As an exam- 
ple, the Sobolev smoothness of a function in the space of bandlimited signals is infinity. 
Algorithms [9j [TUl [12] or even Matlab routines [11] are given for calculating Sobolev 
smoothness of a refinable function. However, those algorithms are not applicable to calcu- 
late the Sobolev smoothness of a function in the space Vp. We give a characterization in 
the next theorem about the Sobolev smoothness of the functions in the space Vp. 

Theorem 5.2. If f G Vp and for some set, 




which is equation (15. 2p . 



□ 



(5.5) 




(5.6) 




then the Sobolev smoothness of the function f satisfies 1/2 (/) > s. 
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Proof. By equations ( II. 2p and ( 13. 2 p we obtain that 



Noting when £G/ n , n<|£|<n+l for each n G Z + , which implies that 1 + £ 2 < 3n 2 , for 
n G Z + . Consequently, we deduce that 

/i/(e)i 2 (i + iei 2 ) s dc < Y.*%r? s I (j2f(^)e- ik A d£ 



where in the last equality again we have used the orthogonality of the set {e lk7T ^ : k G Z} 
on /„, n G Z + . Now apply equation (13.31) we obtain that 

i/(oi 2 (i + \er ^ < 3 i^k^ e„,z + ^ 2s < oo 

11 N ;2 (z+ ) 

by the assumption (15.61) . □ 

For example, consider the space Vp associated with the analytic function F in equation 
( 12.11) such that the sequence b is a geometric sequence with b n = a n (l — a), n G Z + , where 

oo 

a G (—1, 1). The series ^^n 2s a 2 ™(l — a) 2 converges for every s£l. Therefore, we have 

n=0 

the Sobolev smoothness 1*2 (/) = 00. In particular, when a = 0, the space Vp degenerates 
to the space of bandlimited functions, and its Sobolev smoothness is infinity. 

If b is given by b n = ^jjp , n G Z + , the series ri 2s — — - converges if and only if 

n=0 

6 — 2s > 1, that is, s < 2.5, which implies that 1^2 (/) — 2.5 in this case. 
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